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Abstract 

In this paper we initiate the study of SSB in 3 + 1 dimensional rotating, charged, asymptot- 
ically AdS black holes. The theory living on their boundary, Rx S 2 , has the interpretation 
of a 2 + 1 dimensional rotating holographic superconductor. We study the appearance of a 
marginal mode of the condensate as the temperature is decreased. We find that the tran- 
sition temperature depends on the rotation. At temperatures just below T c , the transition 
temperature at zero rotation, there exists a critical value of the rotation, which destroys the 
superconducting order. This behaviour is analogous to the emergence of a critical applied 
magnetic field and we show that the superconductor in fact produces the expected London 
field in the planar limit. 
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1 Introduction 



Over the past year or so, there has been considerable activity in the field of AdS/CMT, which 
attempts to bring a phenomenological approach to the AdS/CFT correspondence to bear 
on problems in condensed- matter physics [IIEIEIIIIIHJEICFJIE]. The scope of this program 
is broad, dealing with such diverse problems as the Nernst effect at a superfluid-insulator 
transition, the quantum Hall effect and cold atoms in the unitarity limit. 

In [9], Gubser suggested that spontaneous £7(1) symmetry breaking by bulk black holes 
- specifically Reissner-Nordstrom-AdS black holes - can be used to construct gravitational 
duals of the transition from normal state to superconducting state in the (not further speci- 
ficed) boundary theory. This can be done by studying the dynamics of a minimally coupled 
complex scalar field with Lagrangian density 



where T>^ = (<9 M — ieA^)^/ and the scalar field mass respects the Breitenlohner-Freedman 
(BF) bound [10] 



In the remainder of this paper, we only consider scalar masses at or above the conformal 
value of m\ = —2. In the first instance one works in the approximation where the scalar does 
not backreact on the geometry. This approximation to the full system is however sufficient 
for an exact determination of the transition temperature of the holographic superconductor. 

Gubser showed numerically that for sufficiently large charge e of the condensate the 
scalar field on the black hole background allows a marginal mode. This is indicative of an 
instability towards a 'hairy' black hole solution with a finite charged condensate outside the 
horizon. The scalar hair breaks the U(l) gauge symmetry in the bulk space time. Hartnoll, 
Herzog and Horowitz [TT] brought this idea to fruition by numerically constructing the state 
with broken symmetry. They were able to show that the electrical resistance in the dual field 
theory indeed drops to zero in the broken phase by studying the fluctuations of the Maxwell 
field around the state of broken symmetry. 

In subsequent works, these authors, and others [TTl[T2"] . verified a number of other physical 
properties of these backgrounds and accumulated evidence that these are gravitationals dual 




(1.1) 



m%L 2 > -9/4 . 



(1.2) 
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of a type II superconductor. This is despite the fact that a local bulk U(l) symmetry 
is a priori associated to a global symmetry of the boundary theory, and hence one would 
expect to find the physics of a charged superfluid [13]. However one may weakly gauge the 
boundary £7(1) and in this regime the material is described by the London theory [11] of 
superconductivity. 

In this paper we aim to put to explore another feature of superconductors, the London 
moment, and by doing so, extend the present framework of holographic superconductors to 
the case of rotating black holes in the bulk. At first sight, this also seems to be a hopeless 
task, since the holographic superconductor does not have a dynamical photon. However, the 
work of [H] showed that it nevertheless produces the screening currents necessary e.g. for 
the Meissner effect. Thus we may also be optimistic about being able to observe the London 
moment. 

At zero angular momentum, the system studied here reduces to a subset of those con- 
sidered in [9], with spherical horizon, where the transition to a superconducting state is 
already known to occur. We thus expect the marginal mode to exist for small values of 
the rotation parameter. An interesting question is whether this mode always persists for 
larger values. Physical intuition tells us that it should not: we know from previous studies 
[12] that a holographic superconductor exhibits a critical magnetic field, above which it is 
energetically favourable to be in the normal phase. We also know that a rotating supercon- 
ductor dynamically generates a magnetic field, the so-called London field p3]. Intuitively 
this arises because the rotation induces a lag between charge carriers close to the surface and 
the charged superfluid in the bulk of the superconductor. Thus we should expect a given 
marginal mode at zero rotation to disappear at some critical rotation. Or, in other words, 
we expect the critical temperature to drop as the angular momentum of the dual black hole 
is raised. We shall demonstrate this behaviour. 

The paper is organised as follows. In section 2 we introduce the dual gravity background, 
the Kerr-Newman-AdS black hole and explain the influence of the bulk rotation on the 
boundary theory. In section 3 we study the angular and radial behaviour of marginal modes 
of a charged complex scalar field in this background and construct the phase diagram of 
the holographic superconductor. For a tachyonic scalar in the planar limit, i.e. close to 
the poles, we compute the London field explicitly. Section 4 is devoted to a discussion and 
appendices A and B contain details on the separation of the wave equation as well as the 
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numerical methods employed in this paper. 

2 Rotation and AdS/CMT 

In this section we introduce the gravitational background, which we argue is dual to a 
two-dimensional material living on the surface of a sphere, exhibiting superconductivity at 
sufficiently low temperature. The study of rotation in the AdS/CFT context was initiated 
in [15] (see also [T6]). 



2.1 AdS4 backgrounds with Angular Momentum 

The desire to introduce rotation into the AdS/CMT story leads us to consider the four- 
dimensional Kerr-Newman-AdS solution, orginally found by Carter [17], who was interested 
in finding space times in which the Hamilton- Jacobi equation would separate. In Boyer- 
Lindquist coordinates, the metric is 



ds 2 = 

P 2 



a • 2,1, /1 2 P 2dr2 P W & d sm 2 8 
d* - - sin 2 Odd) + — K — + F 



with gauge field field 



A = - q f\dt 



A fi 



a sin 2 8 , 
— = — d<p 



and 



adt 



r 2 + a 2 



(2.1) 



(2.2) 



A r := (r 2 + a 2 ) (l + r 2 L- 2 ) - 2Mr + q, 
P 



2 

e > 



J ■— r 2 + a 2 cos 2 8 



A e := I -a 2 L- 2 cos 2 8 
E :=l-a 2 L- 2 (2.3) 



There is some confusion in the literature as to how exactly the parameters a and M are 
related to the physical angular momentum and mass-energy of the black hole, respectively. 
This confusion is cleared up in the work [18] (see also [IH]), by careful thermodynamic 
considerations, and we adopt their definitions 

Ma 



E 



M 

^2 



J 



(2.4) 



for the energy E and angular momentum J of the background. Note that both E and J 
diverge as aL~ x approaches unity. However, J is always strictly bounded above by EL, as 
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expected for a rotating black hole. Finally, we quote the result for the Hawking temperature 
of the black hole [20], 



r+(l + a 2 L- 2 + 3r 2 + L- 2 - (a 2 + q 2 )^ 
47r(a 2 + r 2 



Th = ^— _,:,_ 2 / ^' + ' , (2-5) 



which, in the usual way, is identified with the temperature of the dual field theory living 
on the conformal boundary of the AdS black hole spacetime. The boundary has topology 
R x S 2 , so we are considering a two-dimensional superconductor living on the surface of a 
sphere. The quantity r + is the horizon radius, defined to be the largest real root of the 
equation A r = 0. 

2.2 Rotation and the Boundary Theory 

Massive rotating bodies exhibit the frame- dragging effect: an inertial reference frame outside 
the body is set into rotational motion. This effect diminishes as the distance to the rotating 
body increases. For the original Kerr black hole, which is asymptotically flat, this effect 
vanishes at infinity. However, for the Kerr-AdS family, it does not [18] . Its boundary metric 
is in the conformal class of the Einstein Static Universe 

ds 2 = -dt 2 + L 2 (d9 2 + sin 2 9d4> 2 ^) (2.6) 

where the new angular coordinates satisfy 

4> = 4> — %t , tan# = VE tan 6 . 
L 2 

From the first of these see that the angular momentum of the boundary theory is 

fioo = • (2.7) 

We learn that the local speed of rotation at the equator of the boundary S 2 of radius L reaches 
the speed of light when a = L. In fact the gravitational background becomes singular at 
that point. We therefore restrict to a < L in this paper. 

The meaning of rotation in the field theory becomes apparent if we compute for example 
the free partition function, i.e. the free energy of the theory at zero coupling. It is given by 
the expression 

oo £ 

F = - In Z = T bdry J2 J2 ln (! - m^' m noo) ) , (2-8) 



where i labels the different (particle) species, rji = ±1 for bosons (fermions) and (3 = 1/Tbdry 
The allowed frequencies u follow from the corresponding free wave equations. The summation 
is over the angular momentum quantum numbers of the fields (£, m) as indicated. The 
quantity enters as a chemical potential constraining the angular momentum. 

3 Superconducting Instablity 

We now turn to the main objective of this work, which is to identify an instability under 
perturbations of a charged scalar field propagating on the background. We rely on the 
separation properties of this equation, but relegate the details of the separation procedure 
to appendix A. It is not necessary to follow these in order to understand this paper, but we 
have included them for the interested reader. 

3.1 Radial Equation 

As is well known from a classical result due to Carter [17], the Klein-Gordon equation is 
separable on the background of the Kerr-AdS black hole in 3 + 1 dimensions. This result 
can be extended to the case of a complex (charged) scalar field on the Kerr-Newman-AdS 
background above. We assume that *&(t, r, 8, <fi) ~ e~ luJt ~ m ^ R(r) S (8) , and look for a marginal 
mode, for which m = uj = 0. As in [9|[TT]. it is consistent to seek real solutions R(r) and S(8). 
However, since the complex scalar \I/ is not a gauge-invariant quantity, this is essentially a 
choice of gauge, the details of which are explained in appendix A. For later convenience, let 
us define the horizon function 



where we have set the horizon radius to unityj, and A is the separation constant between 

S(8) and R(r). In the Schwarzschild-AdS limit (a — > 0) it reduces to the familiar value 

£(£+ 1), i G Z being the principal eigenvalue of the spherical harmonics. Sometimes we shall 

lr This can always be done because of scaling symmetry present in the background. However, physically, 
this is nothing but a choice of units in which we measure lengths in multiples of the horizon radius. 



r 2 h{r) = A r (r) . 



(3.1) 



Then we obtain the radial equation 




0. 



(3.2) 
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find it useful to label the values of A by the corresponding integer i to which they reduce. 
While there exists no analytical expression for A, it is possible to determine the numerical 
value of this separation constant by treating the ^-equation as an eigenvalue problem for A. 

To gain a more intuitive understanding of the physics involved in finding the marginal 
mode, let us manipulate expression (I3.2p into the form of a Schrodinger equation. Then 
the question of finding the gravitational mode is re-expressed in terms of finding a marginal 
bound state of the corresponding quantum-mechanical problem. To this end, define the 
tortoise coordinate 

dr 



h(r) 



(3.3) 



Outside the horizon at r = 1, r* is real and takes valuesj r* G (— oo,0). Then the function 
Z(r) = rR(r) satisfies the Schrodinger equation of a zero-energy particle in a potential: 

d 2 Z 



drl 



-V[r m (r),X] Z(n) =0. 



(3.4) 



The potential is given implicitly by 



V(r, A) = h{r) 



2^2 / 



m fl 



eg, 
h(r) 



- - 1 

r 



(h(r))' A + a 2 m% 



(3.5) 



where the prime indicates differentiation with respect to the original radial coordinate r. 
Since we cannot invert r*(r) analytically to get the potential as a function of r, we cannot 
give an explicit functional form of the potential as a function of the tortoise coordinate 
r*. It is however easy to obtain plots of the potential, numerically or otherwise, which are 
sufficient to develop the intuition we wish to achieve. Before proceeding to obtain these plots, 
we must first return to the problem of the separation constant and the angular equation. In 
the process we will gain insight into the nature of the condensate. 



3.2 Angular Equation and Localization 

Without further ado, here is the angular equation for our marginal mode: 
df)(sm9\/Af)d, 



sin 6*vAfi 



+ y/A (dey/A e )do + a 2 sin 2 Qm% + A 



S(6) = 



(3.6) 



2 The upper bound can be seen from the neat identity Xiroots Q'(x ) = ^ ^ or an y polynomial Q(x) at least 

lontou: 

by analytically continuing Q into the complex plane. 



quadratic in x. A simple proof |21j can be given via contour integration of the function f(z) — q^jj defined 
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Again, a detailed derviation is given in appendix A. We can gain some understanding of 
this equation by considering two limits. Firstly, in the non-rotating limit, this equation 
reduces to the associated Legendre equation, defining the 6*-dependent part of the spherical 
harmonics Y£ m (6, <fi). In this limit, we learn that A = £(£ + 1) and it is a useful check on our 
computations to note that our results reduce to these integer values in the appropriate limit. 




Figure 1: Lowest angular AdS spheroidal harmonics. Regions with high density of 
the condensate are shown lighter and regions of low density are shown darker. Both 
panels correspond to a rotation parameter of a = 0.9. On the left panel, we have chosen 
a positive mass-squared term, while on the right the mass-squared term is negative. 

In the flat-space limit, L — > oo, the angular equation reduces to the spheroidal equation, 
first obtained in the black- hole context by Teukolsky [25] . In this limit, no analytical expres- 
sions are known for A, but its values are tabulated. Again, it is a confidence-inspiring check 
on our procedures that our values for A agree with those of the spheroidal harmonics in the 
appropriate limit. For the sake of terminology, let us refer to the solutions of the angular 
equation (13.61) as AdS spheroidal harmonics. The deviation of the AdS spheroidal harmonics 
from being spherical harmonics is measured by the ratio of specific angular momentum of 
the black hole to the AdS length, so we introduce the deformation parameter a — f. 

The eigenvalues of the deformed spheroidal equation (13.61) are then labelled X(£, a, a 2 m\). 
Solutions to equation (I3.6p . much like the associated Legendre functions, have definite parity, 
which allows us to develop a simple, iterative, shooting technique to solve both for the 
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eigenfunction and the eigenvalue to high accuracy. The details are relegated to appendix B. 
There we also list a range of eigenvalues A, corresponding to various values of the deformation 
parameter a. 

An interesting feature of the AdS spheroidal harmonics is that, unlike their French 
cousins, the Legendre functions, the lowest mode is not constant. Rather it has a non- 
trivial profile over the sphere, controlled by the deformation parameter a causing even the 
lowest mode to be localized either near the poles or near the equator. They are localized 
near the poles if m\ < and near the equator if m\ > 0. This is illustrated in figure [U, 
which shows density plots of the lowest mode on the unit S 2 . For the special case of m# = 
they are constant. 

We now have all the ingredients needed to obtain the critical behaviour of the rotating 
superconductor. We start by exhibiting marginal modes of the charged scalar field. 

3.3 Marginal Modes 

3.3.1 Qualitative Considerations 

Evidently, equation (13 .4[) cannot have bound state solutions, unless the potential V(tv) 
develops a negative well in a certain range of parameters. Therefore we will expect the 
lowest marginal mode to lie in the sector with £ = 0. Figure [2|i illustrates this point. The 
occurrence of a negative potential well is closely related to the fact, pointed out by Gubser 
[9], that the effective mass of the scalar field 



gets a contribution due to the coupling to the background gauge potential that can make 
this sufficiently negative for a sufficiently long interval outside the horizon. 

The criterion for instabilty can be illustrated by considering the toy model of the semi- 
infinite square well of depth U and width w. The system will be stable if the potential 
admits negative-energy bound-state solutions and will become unstable once the last of these 
boundstates (i.e. the ground state) acquires positive energy and becomes non-normalizable. 
Elementary quantum mechanics tells us that the condition that there exists at least one 
bound state is 




(3.7) 



where we have briefly reinstated units of h. Thus if the potential well is too shallow or too 
narrow, no bound-state can exist. In particular, for the toy model, a bound state exists for 
v < 4w, where v is the de-Broglie wavelength of a particle of energy U. In searching for 
the zero-energy marginal bound state, we precisely capture the moment where the ground 
state exits the Hilbert space of the theory. We interpret this as the transition point for the 
condensation of the order parameter. 



V(r.) 



V(r) 



1-0 
t-6 








w 

-U 





b) 

Figure 2: a) Plot of effective potential barrier around Kerr-Newman-AdS black hole 
for charged scalar waves with m% = 4. In the chosen coordinates, the horizon is at 
r* = — oo and asymptotic infinity is on the right at r* = 0. b) Square well toy model. 
The depth of the well is U and its width is w. AdS boundary conditions are equivalent 
to putting a reflecting barrier at the origin. 



With this simple picture in mind, let us return to the physics of the holographic super- 
conductor. Keeping the system at at a fixed temperature T below its transition temperature 
for zero rotation, To, the potential well becomes more and more shallow as the specific an- 
gular momentum a is increased. In figure [3] this can be seen as the decrease in depth of 
the potential well, as one increases the rotation a in units of the AdS length. We see that 
the well becomes both less deep and more narrow as the rotation is increased. Hence the 
expectation is that the zero-energy mode will cease to exist at some point. 

This section has illustrated the physical mechanism behind the suppression of supercon- 
ductivity with rotation, but because of the complicated shape of the potential (\3.5L it is 
necessary to analyze the marginal mode in detail, in order to gain a quantitative under- 
standing of the phenomenon. The issue of stability of AdS black holes under scalar field 
perturbations is an interesting subject in its own right. For recent mathematical results on 
(real) scalar fields in AdS black hole backgrounds, see [23] . 
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Figure 3: Plot of effective potential barrier for the A value corresponding to t = and 
eL = 10, m\ = 4. The potential is given as a function of the deformation parameter 
a = aL~ l at constant temperature To. We clearly see that the potential well disappears 
as one increases the rotation. 



3.3.2 Numerical Solutions 

We will now embark on a study of the bevhaviour of the marginal mode as a function of the 
parameter a. It should come as no surprise to the reader that the only successful quantitative 
approach approach to solving the radial equation is the numerical one. Already in the much 
simpler case of RNAdS black holes, no analytical solutions are known for the marginal mode, 
much less for the broken phase. 

In order to facilitate the numerical analysis, let us transform to the variable z = r~ 1 . 
Then the question of the marginal mode is transformed into a boundary value problem on 
the interval (0, 1). We wish to solve the equation 

~ e 2 q 2 e (z-l) 2 _ m%(a 2 z 2 + l) _ V 
z 4 h(z) z 4 z\ 



d / . .dR\ 



R(z) = (3.8) 



subject to certain boundary conditions. We see that both the horizon at z = 1 and asymp- 
totic infinity at z = are regular singular points of this equation. Near the AdS boundary 
at z — 0, the solutions take the familiar form 

R{z) ~ + ^ 2 z A - , (3.9) 

where A± is the conformal weight of the dual operator, given by the root of the indicial 
equation 

A(A - 3) = m\L 2 . (3.10) 
10 



Unitarity requires that A > |. Depending on the bulk scalar mass, either or both of the 
operators 

(Ox) = v^i 

(0 2 ) = V2ty 2 (3.11) 

may condense. We have chosen the same normalization for the bulk-boundary coupling as 
[TT] . The indicial equation at the horizon tells us that the solution there behaves as 

R(z) ~ R {0) + R {1) \n(z - 1) (3.12) 

and regularity demands that R^ = 0. 

In order to solve ( 13. 81) . we set either of = 0, corresponding to either (0\) or ((9 2 ) 
condensing. Thus we develop the solution in a Frobenius series around infinity 

R(z) = z A ^ + a n z n ^j i = 1 or 2 (3.13) 

This can be done to very high order using a symbolic algebra package to solve recursively 
for the coefficients. Similarly the expansion around the horizon takes the form 

i2(z) = l + ( 3 - 14 ) 

n 

i.e. we have chosen to normalize R^ = 1. As remarked in |9| it would be more proper to 
treat R^ as a small parameter that allows us to justify neglecting the backreaction, but in 
the end this is a mere rescaling of the mode. 

We use numerical integration to match the two series solutions at some intermediate 
point by adjusting the coefficient This shooting technique can be implemented^ very 
efficiently using a numerical root finder and looping over initial conditions. 

However, after matching the value of the function, there is no guarantee that the first 
derivatives will also match at this point. This is because there is only a marginal mode at 
certain specific values of the background parameters. Evaluating Th, given in equation ( 12. 51) . 
at the critical values for which the mode first appears gives us the critical temperature. An 
illustration of several marginal modes is given in figure HI 



3 As with all numerical procedures described in this paper, I shall make Mathematica™notebooks avail- 
able upon request. 
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Figure 4: Examples of lowest (i.e. no nodes) marginal modes for condensation of 
0\. We take eL = 8, m 2 L 2 = 4 and a = 0,0.7,0.9 for the solid, dashed and dot- 
dashed modes respectively. In the Schrodinger representation these are the radial 
wavefunctions of the zero-energy bound states. 

The family of backgrounds we are interested in depends on the parameters L, a, q e , r + 
and the scalar field has mass mq, and charge e. We have chosen units such that r + = 1 
and analyse the equation for fixed values of e and m^. This results in a curve of critical 
points in the (a, g e )-plane. For concreteness we fix q e = 1 to extract a value for T c . Thus we 
are looking at projections of the full phase diagram. It would be interesting to extend this 
analysis to explore the entire phase diagram. 

3.4 Phase Diagram 

In order to find the critical temperature at a given value of background parameters, we follow 
the example of [22] and focus on the condensation of 0\ for different parameter values. This 
is convenient because it always corresponds to a normalizable solution, as long as m% > — |. 
Quantitatively, this only gives a lower bound on the critical temperature, since in the mass 
range where both 0\ and O2 are normalizable, 0\ may in fact condense first. However, in 
this paper we are not so much concerned with the actual value of T c , but rather its behaviour 
function of aL . 

We find that in all cases T c decreases as the angular momentum of the superconductor is 
increased. This means that at a temperature slightly below T , the critical temperature at 
zero rotation, there always exists a critical value of aL -1 above which the material is forced 
back into its normal phase. 

However, for some choices of parameters, such as in figure [5^, one can choose the tem- 
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Figure 5: Critical temperature for the condensation of 0\ at q e = 1 as a function of 
a = a/L. Left panel for m\L 2 = —2 and right panel for m|,L 2 = 4. Superconductivity 
is increasingly suppressed as the rotation is increased. 

perature low enough, that no amount of rotation will force the material back into its normal 
state. Recall that for aL^ 1 = 1 the rotation speed at the equator equals the speed of light. 
However, as we have seen in figure [U the condensate may be localised away from the equator, 
close to the poles, where the local speed of rotation is arbitrarily small. Thus, it is plausible, 
that in these cases, enough of the condensate recedes far enough away from the areas of the 
S 2 , where the rotation would destroy the superconducting order. 

3.5 Planar Limit and London Moment 

Let us consider what happens if we concentrate on a small region around the north pole of 
the boundary S 2 . To this end, we introduce the coordinate 



and consider the regime u <C L, with a — t held constant. In other words we zoom in 
on the physics close to the north pole on scales much smaller than the AdS length. In this 
regime the curvature of the sphere is negligible and we are effectively dealing with a flat 
horizon rotating with angular velocity f2 = jk. Note that the local speed of rotation is of 
order O (mL _1 ), small compared to the speed of light. The angular equation (13.61) simplifies 
to give 



u = L6, 



(3.15) 




(3.16) 



13 



Compare this with the equation of a planar holographic superconductor immersed in a 
magnetic field, viz. 

S"(u) + ^S'{u) - (^Pj S{u) = -XS{u) , (3.17) 

for a different separation constant of mass dimension two, as defined in [llj. This equation 
was studied in pj], [12], where it was shown to lead to superconducting droplets that are 
exponentially confined near the origin (in our case the north pole). This behaviour is con- 
sistent with our findings for operators that are dual to tachyonic bulk scalars. We conclude 
that in these cases, the rotating superconductor is equivalent to a static superconductor, 
immersed in a magnetic field of magnitude 




(3.18) 



where £1^ is the angular momentum of the boundary theory and m = \m^,\. Again, this field 
is weak, of order O (£ _1 ), in the limit we are considering. Furthermore, in the nonrelativistic 
regime a 2 L~ 2 <C 1, it reduces to the relation between angular momentum and magnetic field 
first obtained by F. London [TJ], by considering his phenomenlogical theory for a rotating 
superconductor. The work of [25] studied the case where one adds linear momentum to the 
system. It would be interesting to see if their results can be seen in the present framework 
by zooming in on the equator of the solutions. 

4 Discussion 

In this paper we have demonstrated that finite-temperature gravitational backgrounds with 
non-zero angular momentum exhibit SSB much like their static limits. Furthermore we have 
seen that the additional dependence on rotation causes phenomena akin to those found in 
studies of superconductors in magnetic fields. 

An important difference to this case is that there always exists a critical magnetic field 
B c , even at very low temperature, above which the superconducting order is destroyed. We 
have seen that this is not always the case for the rotating superconductor on a sphere. This 
can be explained by noting that the local speed of rotation decreases as one moves away 
from the equator, and thus there can be regions of material, where the superconducting 
order is present, even when the equator moves at the speed of light. Near the poles we 
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have seen that the field produced is precisely that predicted by F. London on the basis of 
his phenomenological theory of superconductors. The field is parametrically weak, of order 
C(L _1 ), so at low temperatures, it is not strong enough to destroy the superconducting 
order. This fact is also borne out in the full numerics, as exemplified in figure 5a. 

The London field is often written in terms of the dressed mass m* of the condensate, 
rather than the bare mass m. In our case, we find that m* = -^=. Thus the dressed mass 
corresponds to the relativistic mass of a particle of mass m at the equator of the boundary 
sphere. 

We have seen that even the lowest angular eigenfunction localizes the condensate in a 
droplet close to the equator for positive mass squared of the dual scalar field, and in a ring 
around the equator for negative mass squared. It is interesting that the localization behaviour 
of the condensate is qualitatively different for the case of positive bulk mass compared to 
negative mass (squared) and only in the latter case have we offered an explanation in terms 
of the material's London moment. It would clearly be an interesting question to understand 
what the different localization properties mean from the point of view of the boundary 
superconductor. It would appear that in one case the instablitity is towards forming a 
vortex anti-vortex pair localized on the antipodes, while in the other case we have a pair of 
lumps of superconducting material on the two poles. 

Our analysis was done to first order in the small condensate. The localization behaviour 
observed in this paper would be expected in the early stages of vortex formation. At the 
level of approximation of this paper, the model does not allow for dynamically generated 
screening currents. However, (TTJ [12] argue that going to higher order in the condensate will 
generate the screening currents necessary to confine the normal phase in vortex cores whose 
centres will allow magnetic flux to penetrate. Thus our results confirm the expectation that 
the rotating holographic superconductor has a London moment generating a London field. 

It would be very interesting to extend the computations presented in this work to include 
the effects of backreaction of the scalar, both on the Maxwell field and on the gravitational 
background in order to understand the formation of vortices from the localized droplets in 
our system and to directly see the London field. 

Independently of any holographic interpretation, the results of this paper indicate the 
existence of a new branch of stationary black hole solutions. Here we have found a marginal 
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mode of the stationary asymptotically AdS black hole, which preserves axial symmetry. 
Just like in the static case, where the existence of a static marginal mode preserving the full 
spherical symmetry indicates the existence of a new branch of charged static solutions, here 
we expect the existence of a new branch of charged stationary solutions. Clearly these last 
two issues are numerically more involved and will require more sophisticated methods than 
the ones discussed in this paper. 
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and comments. I gratefully acknowledge generous support by EPSRC and Trinity College, 
Cambridge. The work presented here was initiated at MIT while I was supported by the NSF 
under grant PHY-0600465 and the DOE under collaborative research agreement DEFG02- 
05ER41360. 

A Separation of Variables 

While it can be shown by direct computation that the charged scalar field allows separable 
solutions, this procedure is as tedious as it is inelegant. We shall instead use the dyadic 
index formulation of Newman and Penrose [26] which is well suited to the present problem. 
A suitable Newman-Penrose tetrad is via the differential operators 



Acknowledgements 



D = Fd t 



((r 2 + a 2 )d t + A r d r + aZdtf,) 



■r 



A = n^d, 



- ((r 2 + a 2 )d t - A r d r + aZd^) 




(A.l) 



where all quantities are as defined in the main paper and in addition 



p = r + ia cos 6 



(A.2) 



As ever, the metric takes the form 



ds 2 = 2{m®m — £<g)n) 



(A.3) 
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in terms of the one-forms dual to (1A.1I) . A simple representation of the non- vanishing spin 
coefficients is given by 



7i = — <5 hip* 
p = D\np* 

H = -A In p* 

1- f asm9^/AgA 1 

a = -0 In ; tt; 

2 V (P*) 2 



(A.4) 



r = 5 In p* , 

13 = --Shy (a 2 sm 2 6A e ) 

1 A1 ( asin^v^AeA, 
7 = - A In ; — — 



(A.5) 



where, in keeping with the notation introduced by Newman and Penrose, we also use the 
letter p for one of the complex spin coefficients. Finally, the charged scalar field obeys the 
equation 



V 



(a)(6) 



(V (a) + ieA (o) )(V(6) + zeA (b) )^ = 



(A.6) 



where we use Chandrasekhar's notation of putting indices that refer to the NP tetrad in 
brackets. Let us introduce the family of differential operators 

iK 

T> n = d r — - — h nd r In A r 

iK 

T>1 = d r + -: h nd r In A r 

A r 



Cl = do 



H 

a7 

H 



+ ndg In sin B*J Ag 



£ n = do — — — h nde In sin 9y Ag 
Ag 

where acting on functions of the form 

ty(t,r,6,<f)) = e- l(uJt+m ^y(r,9), 

the constants K and H are 



K = (r 2 + a 2 )c<j + amS , 



H 
AT 



wa sm 6/ H 

sm # 



(A.7) 



(A.J 



(A.9) 



It is then straightforward to show in a few steps that the charged scalar equation simplifies 
to the expression 

ieq e r 



r a, 

w 



A, 



A r 



A r 



A, 



£ (4 Va^a, + £iVa74) } $ = m|$ 



(A.10) 
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Now, notice that the operators X>o,i an d their hermitian conjugates are purely radial, while 
the operators £0,1 and their hermitian conjugates are purely angular. Thus the equation is 
separable upon multiplying by p 2 , resulting in the two equations 



& r (V{Dl + V\V ) + 2ieq e r(V - V\) + 2 



e 2 q 2 r 2 



2{r 2 + a 2 )m% - 2A 



R(r) = (A.ll) 



and 



(v^£iv^4 + V^4v^A)) +2a 2 sin 2 #m| + 2A 0(9) = 0. (A. 12) 

Substituting the definitions of the differential operators £0,1 and Po,i results in the radial 
and angular equations quoted in the main text modulo a gauge choice, on which we now 
elaborate. 



B Details on Numerical Methods 

This appendix contains details on the gauge choice used in equation (13. 2p and gives details 
on the numerical algorithms used to compute the separation constants used in the bulk of 
the paper. 

B.l Gauge Choice and Regular Solutions 

The potential contribution to the radial equation contains an effective mass-squared term 
proportional to the square of the gauge field. This term seemingly blows up at the horizon 
unless the gauge field A goes to zero there. Notice however, that the function R(r) that we 
are solving for is itself not a gauge invariant quantity. In particular, it has a gauge-dependent 
phase. Thus if we choose the radial function to be 

e ™V{r)R(r) , (B.l) 

we can remove the apparent singular behaviour of the potential at the horizon with a par- 
ticular gauge choice for <p(r). A computation shows that the function 

¥>(0 = — / TTT = > ( B -2) 

r+ J h[r) r + 

which we recognize as the tortoise coordinate f 1 3 . 3 j) . removes the apparent singular behaviour 
at the horizon. With this choice of phase, the gauge-invariant part of the radial function 
satisfies the equation (13. 2\ . 
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B.2 Angular Eigenvalues and Eigenfunctions 

The separated equations are coupled through the separation constant A. It is thus important 
to have precise numerical methods to determine it. With x = cos 9 equation (13. 6p reads 



a 2 x 2 



d_ 

dx 



\/l — a 2 x 2 (l — x 



2\ ULJ lm 



dS. 



+ 



dx 

A + a 2 m\{l — x 2 ) ■ 



dS, 



- a 2 x(l - x 2 )^^ 
dx 

m 2 (l — a 2 



x 2 



Sj m {x) = 0, (B.3) 



where we have briefly reinstated the azimuthal quantum number m. In the limit L — > oo, 
equation (1B.3j) goes over into the spheroidal equation and the eigenfunctions S] m (x) are the 
spheroidal harmonics. In many references on spheroidal harmonics the parameter a 2 m|, is 
called 7 2 . In the limit a — > it becomes the associated Legendre equation. 

In this paper, we only consider Sj m (x, A) for m = 0, which is all we need in determining 
the critical temperature. Note, however, that the angular eigenvalue of the (AdS) spheroidal 
equation depends on m, unlike the more well-known case of the spherical harmonics (i.e. the 
Legendre polynomials). However, our procedures work equally well for non-zero values of m. 

B.2.1 Numerics 

We now describe the numerical procedure used to extract the eigenvalues A from this equa- 
tion. It is evident that solutions of equation ( IB. 31) are either even or odd functions. We will 
make use of this crucial fact, which makes it easier to use a simple shooting method: 

1. Obtain series expansions of Sj m at x = ±1. There is a regular branch and a logarithmic 
branch of solutions. As usual we select the regular branch. We set the following 
boundary conditions 

S lm ( - 1 ) = S lm ( 1 ) = 1 S Jm ( X ) is eVen 

52n(-l) = -SLV) = 1 Sj m {x) is odd (B.4) 

The choice of the magnitude of the eigenfunctions at either end of the interval is clearly 
just a choice of normalization. 

2. Numerically integrate inwards to the origin from both ends for a fine grid of values 
of the eigenvalue A in order to match the Frobenius expansions at either end of the 



19 



intervals. This gives rise to a series of functions S] m (x, A), where the ± refers to the 
numerical integration from ±1. 

3. For odd functions the first derivatives of Sj m (x, A) are guaranteed to coincide at the 
origin. Thus we define the function 

O(X) = S]^(0,X)-S^(0,X) 

by interpolating the values of the discrete grid of values of A. For even eigenfunctions, 
we are guaranteed that the value of S] m (x, A) matches at the origin. Thus we define 

E(X) = (S^)'(0,X)-(S]-)'(0,X) 

4. Find the zeroes of -E'(A) and 0(A). Every root corresponds to an eigenvalue of the 
angular equation. 

We now tabulate (a subset of) the A values that were used to find the critical temperatures 
displayed in figure [51 Table 1 lists a representative sample of eigenvalues to five significant 
digits that were used in obtaining the phase curves of critical temperature versus rotation. 
Our numerical algorithms in fact allow a far higher accuracy, but it would be impractical 
and of little value to list more significant digits here. 

Table 1: Angular eigenvalues of AdS spheroidal harmonics. The first line gives the 
deformation parameter a = aL~ l and the following two lines supply the angular eigen- 
values for the two mass-squared values used to compute the phase diagram. 



a 





0.1 


0.2 


0.4 


0.6 


0.8 


0.9 


m%L 2 = -2 





0.0133 


0.0532 


0.2115 


0.4690 


0.8086 


0.9921 


m%L 2 = 4 





-0.0266 


-0.1071 


-0.4340 


-0.9997 


-1.8393 


-2.3764 
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